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We determine the distribution of linearly polarized gluons of a dense target at small x by solving 
the B-JIMWLK rapidity evolution equations. From these solutions we estimate the amplitude of 
~ cos 20 azimuthal asymmetries in DIS dijet production at high energies. We find sizeable long- 
range in rapidity azimuthal asymmetries with a magnitude in the range of V 2 = (cos 20} ~ 10%. 


Transverse momentum dependent (TMD) factoriza¬ 
tion [1, 2] in deep inelastic scattering predicts a distribu¬ 
tion for linearly polarized gluons in an unpolarized tar¬ 
get [3, 4]. This is reflected in cos 20 asymmetries in dijet 
production [5, 6] and in other processes [7-9]. To date 
little is known about the magnitude of these functions in 
the small-# regime of high energies. In this paper we per¬ 
form first estimates of these functions by solving the B- 
JIMWLK renormalization group equations [10-21]. Also, 
we use our solutions to analyze the magnitude of the re¬ 
sulting cos 20 asymmetry in dijet production [5, 22] at 
leading order. These could be tested at a future electron- 
ion collider (EIC) [23, 24], where the small-# effects dis¬ 
cussed here can be enhanced by using a nuclear target. 

Recent data for high multiplicity p+p [25, 26] and 


p+Pb [27-33] data at the LHC have revealed long-range 
(in rapidity) angular cos 20 “ridge” correlations in par¬ 
ticle production high multiplicity events. The magni¬ 
tude of these long range correlations is conventionally 
parametrized in terms of #2 = (cos20). In fact, the 
azimuthal correlation in DIS dijet production at high 
energy originates also from the long-ranged eikonal in¬ 
teraction and so results in a similar experimental sig¬ 
nature as the “ridge”. To make this connection ex¬ 
plicit we shall parametrize the azimuthal structure aris¬ 
ing from the linearly polarized gluon distribution in terms 
of #2 = (cos 20), and determine its dependence on the ra¬ 
pidity imbalance of the dijet. 

At leading order the cross section for inclusive produc¬ 
tion of a dijet in 7 *-nucleus scattering is given by [5, 6] 
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Here, 

P± = (1 - z)ki - zk 2 , q± = ki + k 2 (3) 
are the dijet transverse momentum scale Pj_ and the 


transverse momentum imbalance q±, respectively. The 
transverse momenta of the produced quark and anti¬ 
quark are given by ki and k 2 and their respective light- 
cone momentum fractions are z and 1 — z; the dijet in- 
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FIG. 1: Linearly polarized and unpolarized WW gluon distributions versus transverse momentum q± at different rapidities Y. 
Transverse momentum is measured in units of the saturation momentum Q s (y). The curves correspond to evolution at fixed 
a s = 0.15. 


variant mass is given by M = P_l/ sj z(l — z). Also, 
e 2 = z( 1 — z)Q 2 with Q 2 of order P]_. Here, we re¬ 
strict ourselves to kinematic configurations where P± is 
greater than gj_, referred to as the “correlation limit” in 
Refs. [5, 22], 

In Eq. (2) </> denotes the azimuthal angle between P± 
and q i , respectively. We introduce the following measure 
for the azimuthal anisotropy, 

v 2 = (cos 2 <j>) . (4) 


The average over (f) in this equation is performed with 
the weight (1) or (2), respectively. Since 
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( 5 ) 


is independent of <f >, for a longitudinally/transversally po¬ 
larized photon we have 
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( 6 ) 

The linearly polarized and unpolarized dis¬ 
tributions are dehned as the traceless part and the trace 
of the Weizsacker-Williams unintegrated gluon distribu¬ 
tion, respectively: 


and cross sections for different observables can be related 
to different correlation functions of the Wilson lines. The 
Wilson line is a path ordered exponential of the covariant 
gauge field, whose largest component is A + : 

UM-T^igJte-A^-.xr)}. ( 8 ) 

The Weizsacker-Williams unintegrated gluon distribu¬ 
tion [5, 22, 34], on the other hand, is expressed most 
naturally in terms of the light cone gauge ( A + = 0) field, 
which has large transverse components. These can be 
obtained by a gauge transformation 

A*(x t ) = — I7 t (x T )5 i 17(x T ) . (9) 

Since, in light cone gauge, the gauge held lives above 
the light cone A 1 (xt,x~) ~ 0(o; _ )A i (xt), this held can 
also be thought of as a sheet of color electric held on the 
light cone E 1 (xt, x~) = S(t— z)A*(xt). The Weizsacker- 
Williams distribution is simply the two-point correlator 
of the light cone gauge helds 
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xGl w = l&xGW - I - 2^ xhf . (7) 

In the CGC framework the gluonic degrees of freedom at 
small x are described by Wilson lines. They are path or¬ 
dered exponentials in the strong color held of the target, 


where we have normalized the distribution with the 
transverse area of the target L 2 . This normalization 
drops out of the results expressed in terms of the ellip¬ 
tical asymmetry v 2 . For analytical calculations of the 
functions <7_l) and h^\xo,q±) in the McLerran- 

Venugopalan (MV) model [35, 36], see Refs. [6, 22]. 
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FIG. 2: The average azimuthal anisotropy V 2 = (cos 20} versus the dijet transverse momentum scale Pt or the dijet transverse 
momentum imbalance qr, respectively. The assumed 7 * A center of mass energy is y/s = 100 GeV. Since Q 2 = 4P 2 and z — 1/2 
these curves apply to either longitudinal or transverse photon polarization. Solid (dashed) lines correspond to fixed (running) 
coupling evolution. 


We obtain the Wilson lines U numerically from B- 
JIMWLK evolution in Y = \n.(xo/x), starting from 
an initial condition at xq = 10 ~ 2 using the the MV 
model. The initial condition on the lattice is constructed 
as described in detail in Ref. [37]. The B-JIMWLK 
equation can be solved on the lattice with a Langevin 
method [38, 39]. We use here the “left-right” symmet¬ 
ric [40] numerical method introduced in Ref. [41], using 
either fixed coupling or a running coupling with the al¬ 
gorithm of Ref. [41]. As in e.g. Ref. [42], we determine 
the saturation scale Q s numerically from the two-point 
(dipole) function of the Wilson lines. The renormaliza¬ 
tion group evolution increases Q s roughly as Q 2 ~ a; -0 ’ 3 . 
For the calculation of the light cone gauge field one needs 
Fourier transforms of derivatives of Wilson lines. Some 
care must be exercised to obtain the proper momen¬ 
tum space distribution: we have used two different cen¬ 
tered difference methods (discretizing over one or two 
lattice spacings) and found that the results are equiva¬ 
lent. For the fixed coupling evolution we take a s = 0.15 
to provide an evolution speed roughly in line with in¬ 
clusive HERA data. For running coupling we use in 
this preliminary study the slightly overestimated value 
Qs( 2 ;o)/Aqcd = 11 , which also slows down the evolution 
closer to experimentally observed values. 

For our numerical estimates below we take Q 2 = 4P]b 
Hence, for z = 1/2, v% and v have equal magnitude 
but there is a relative phase shift of n/2. The physi¬ 
cal momentum scale is set by the saturation momentum 
at Xq. To obtain the numerical values in the plots we 
take Q s (xo) = 1 GeV (for a qq dipole). The saturation 
momentum corresponds to the scale where the forward 


scattering amplitude is of order 1 . 

We now turn to describe our results. We first show the 
solution for the unintegrated gluon distributions before 
discussing the azimuthal asymmetry w.r.t. the direction 
of q _1 of the 7 * A cross section. 

Figure 1 shows the dependence of G*- 1 -* and at dif¬ 
ferent evolution rapidities Y on transverse momentum. 
We refrain from showing curves for running coupling evo¬ 
lution since they look very similar. Either one of the 
TMDs drops rapidly as a power of q± at high transverse 
momentum q± Q s and so they are best measured at 
q |_ of order a few times the saturation scale. For a heavy- 
ion target the saturation scale is boosted (on average over 
impact parameters) by a factor of ~ A 1 / 3 [43] which fa¬ 
cilitates such measurements in a regime of semi-hard q±. 

The degree of gluon linear polarization is maximal at 
high transverse momentum, /G^ —> 1 ; the satura¬ 
tion of the positivity bound of the cross section has also 
been observed in perturbative twist -2 calculations of the 
small -x field of a fast quark [4, 6 ]. On the other hand 
/G^ -C 1 at low q± which conforms to the expected 
power suppression. At fixed q±/Q s (x) the ratio of these 
functions decreases rather slowly with rapidity, at least 
after an initial evolution away from the MV model to¬ 
wards the B-JIMWLK fixed point. This means that, be¬ 
cause of the growth of Q s , the ratio /G^ at fixed 
transverse momentum q± decreases with rapidity. Thus 
the emission of additional small -x gluons reduces the de¬ 
gree of polarization. Our results show that this effect 
can quite well be parametrized by geometric scaling as a 
universal function of q±/Q s . 

In Fig. 2 we show the elliptic asymmetry as a func- 
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FIG. 3: The average azimuthal anisotropy V 2 = (cos 20} ver¬ 
sus the dijet rapidity imbalance £ = log (1 — z)/z. Thick 
(thin) lines correspond to longitudinal (transverse) photon 
polarization. 

tion of the dijet transverse momentum scale P± and the 
transverse momentum asymmetry q±. Increasing P± in¬ 
creases x and suppresses evolution effects and so V 2 {P. l) 
increases towards the MV model initial condition. The 
reason for the difference between the fixed and running 
coupling curves in V 2 (P±) is that in this preliminary 
study they have not been adjusted to have the same evo¬ 
lution speed dy ln<3 2 (Y). We observe the same behavior 
for even though x increases only slowly with q±; 

here the increase of the elliptic asymmetry is mainly due 
to h ( l ) (q±)/G( 1 '>(q±) -> 1 as q±/Q s 1, as shown above. 
Overall, in the kinematic range considered in Fig. 2 we 
find a rather substantial magnitude of V 2 ~ 10 %. 

Figure 3 shows V 2 versus the rapidity asymmetry 

£ = log — ■ (ii) 

z 

Our calculation applies for moderately large rapidity sep¬ 
arations less than l/a s , since we are assuming that the 
two jets are sensitive to the same distribution of Wilson 
lines. We find a mild increase of V 2 away from z = 1/2 
which is due to the fact that asymmetric dijet configura¬ 
tions probe the gluon field of the target at larger values 
of x. The slow evolution of the eikonal interaction with 
x translates into a rather flat d 2 (£) over several units in 
£ away from the boundary of phase space. Hence, at 
high energies the azimuthal asymmetry is long range in 
rapidity. 

In summary, we have computed the TMD distribu¬ 
tion of linearly polarized gluons for a large nucleus 
at small x. We have used the McLerran-Venugopalan 
model to obtain initial conditions at xq ~ 10~ 2 and 


the B-JIMWLK equations to evolve to lower x. We 
find that for realistic values of x and transverse mo¬ 
mentum imbalance q± that h±\x,q±) is of substantial 
magnitude. This results in large elliptic azimuthal asym¬ 
metries V 2 = (cos 20) ~ 10% in DIS dijet production. 
Also, the azimuthal correlations are long range in rapid¬ 
ity, i.e. V 2 depends weakly on the rapidity asymmetry 
£ = log(l - z)/z. 

In the future we intend to check other initial condi¬ 
tions for the evolution, although we do not expect quali¬ 
tative modifications of the results presented here. It will 
be interesting, also, to study Sudakov resummation ef¬ 
fects [44, 45] as well as more general kinematic configu¬ 
rations which require quadrupole matrix elements [5]. 
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